Quantum-enhanced optical phase tracking 

Hidehiro Yonezawa, 1 Daisuke Nakane, 1 Trevor A. Wheatley, 1 ' 2,3 
Kohjiro Iwasawa, 1 Shuntaro Takeda, 1 Hajime Arao, 1 Kentaro Ohki, 4 

Koji Tsumura, 5 Dominic W. Berry, 6 7 Timothy C. Ralph, 2 8 
Howard M. Wiseman, 2 ' 9 * Elanor H. Huntington, 2 3 Akira Furusawa 1 * 

department of Applied Physics, School of Engineering, 
^ The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 1 13-8656, Japan 

Centre for Quantum Computation and Communication Technology, Australian Research Council 
3 School of Engineering and Information Technology, University College, 
The University of New South Wales, Canberra ACT 2600, Australia 
4 Department of Applied Mathematics and Physics, Graduate School of Informatics, 
Kyoto University, Yoshida-Honmachi, Sakyo-ku, Kyoto 606-8501, Japan 
5 Department of Information Physics and Computing, 
The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan 
institute for Quantum Computing, University of Waterloo, Ontario N2L 3G1, Canada 
7 Department of Physics and Astronomy, Macquarie University, NSW 2109, Australia 
8 School of Mathematics and Physics, University of Queensland, Brisbane QLD 4072, Australia 
9 Centre for Quantum Dynamics, Griffith University, Brisbane QLD 4111, Australia 

*To whom correspondence should be addressed. E-mail: 
h.wiseman@griffith.edu.au (H.W.); akiraf@ap.t.u-tokyo. ac.jp (A.F.) 



Tracking a randomly varying optical phase is a key task in metrology, with 
applications in optical communication. The best precision for optical phase 
tracking has till now been limited by the quantum vacuum fluctuations of co- 
herent light. Here we surpass this coherent-state limit by using a continuous- 
wave beam in a phase-squeezed quantum state. Unlike in previous squeezing- 
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enhanced metrology, restricted to phases with very small variation, the best 
tracking precision (for a fixed light intensity) is achieved for a finite degree 
of squeezing, due to Heisenberg's uncertainty principle. By optimizing the 
squeezing we track the phase with a mean square error 15±4% below the 
coherent-state limit. 

There are many tasks where precise optical phase estimation is critical, including commu- 
nication (1, 2), and metrology (3). Quantum mechanics imposes a fundamental bound on preci- 
sion (4-6), and this already limits gravitational wave detection (7-9) and can guarantee security 
in quantum cryptography (10). The quantum limits are determined by optimizing (subject to 
constraints) the input quantum state, the quantum measurement, and the data processing. Much 
effort has been devoted to approaching the fundamental quantum limits (3, 5, 6). 

Phase estimation can be divided into two kinds (11): phase sensing, where the phase is 
known to always lie within some small interval (e.g. (12)), and general phase estimation, where 
it is not so constrained (e.g. (11)). In the former case, when (as in most practical situations) 
the field has a large coherent amplitude, the problem can be linearized in terms of the phase 
rotation (7), which greatly simplifies the task of optimizing the input state and the measurement. 
By contrast, in the case of unconstrained phase estimation the problem cannot be linearized, 
and as a consequence the optimization problem is considerably harder (11, 13-22). While a 
quantum enhancement of phase sensing using nonclassical states of light has recently been 
demonstrated (8, 9), this has been done for general phase estimation only with post-selected 
results (11). 

We present a demonstration of unconstrained phase estimation with a quantum enhance- 
ment using nonclassical (squeezed) states. We use homodyne detection, with no post-selection 
of data, and no compensation for losses or detector inefficiency in the system. Moreover, the 
problem of a stochastically varying phase is addressed (8, 9, 19-22), as is highly relevant for 
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physical metrology and communication, rather than a time-invariant (but initially unknown) 
phase (11-13, 15-18). To perform optimal estimation we have implemented optical phase track- 
ing — a phase-lock loop which strives to maintain the maximum measurement sensitivity for a 
widely varying phase. The quantum noise in the photocurrent prevents the maximum sensitivity 
from being perfectly maintained, which is why we observe an optimal degree of squeezing. 

Our experiment (Fig. 1A) uses a continuous-wave optical phase-squeezed beam. The phase 
of the beam is modulated with the signal ip(t), the waveform to be estimated (22-24). As in 
Refs. (21, 22), we use a stochastic waveform defined by 



Here dV(s) is a classically generated Wiener process (25) (white noise), A -1 is the correlation 
time of ip(t), and n determines the magnitude of the phase variation, which is of order unity. 
This ip(t) is a continuous-time random walk with a tendency to return to the mean phase of 
zero, a kind of noisy relaxation process that occurs in many physical situations (25). 

The phase-modulated beam is measured by homodyne detection, using a local oscillator 
(LO), yielding a noisy current I{t). The LO phase $(£) is feedback-controlled to be $(£) « 
7r/2 + (p(t), as this is the most sensitive operating point for sensing changes in ip{t) (Fig. IB). 
Because ip(t) is unknown, the best strategy is adaptive metrology (11, 13, 15-21), in which 
feedback control is used to set $(£) = if fit) + n/2, where iff(t) is a filtered estimate of ip(t) 
— that is, an estimate based on I(s) for all s < t. This gives a normalized homodyne output 
current I(t) of (19,26), 



Here |a| is the amplitude of the input phase-squeezed beam, and dW(t) is another Wiener pro- 
cess (25), arising from the squeezed vacuum fluctuations. The magnitude -R sq of these quantum 
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fluctuations is determined by the degree of squeezing (r m > 0) and anti- squeezing (r p > r m ), 
and by aj = {[<p(t) — </?/(t)] 2 ). Note that several approximations — justified in Ref. (26) — 
have been made to derive Eq.|2l most notably, a second-order expansion for I(t)dt in the small 
variable [<p(t) — <ff(t)]. 

For optimal feedback control, the Kalman filter is used for (pf(t) (22), which is the causal 
(i.e. real time) estimator with the lowest mean square error (MSE). The Kalman filter is the 
optimal filter for estimating ip(t) of the form of Eq. \T\ when using a coherent beam (22), and the 
calculation generalizes to our squeezed case (26). Though the filtered estimate ip/(t) is a good 
estimate of the signal phase (f(t), to obtain the best estimate we apply the acausal technique of 
smoothing (21, 22, 24). After storing data over a certain period of time, a precise estimate of 
<p s (t) is obtained at a time t in the middle of that period by using observations both before and 
after t. The MSE of the smoothed estimate a 2 s = (\tp(t) — y2 s (t)] 2 ) is given as (22, 26), 

a 2 s = k/ (2 \]^\a\ 2 /i? sq + A^ . (3) 

Recall that R sq (Eq.O is a function of aj, so the above expression for <r 2 is still implicit. The full 
solutions are given in (26), but in the parameter regime of our experiment, cr 2 is approximately 
proportional to >J R sq . That is, by using a nonclassical beam with effective squeezing _R sq < 1 
we expect to be able to overcome the coherent- state limit (CSL) by a factor of \J l\ q . 

Our experiment (see Fig. 1C) uses an 860 nm continuous-wave Titanium Sapphire laser. The 
phase-squeezed beam is added by an optical parametric oscillator (OPO). The OPO is driven 
below threshold by a 430 nm pump beam, generated by a second-harmonic-generation cavity. 
We obtain up to —4 dB of phase squeezing. The signal (p(t) is produced by a digital random 
signal generator and a low-pass filter with a cutoff frequency of \/2ir. This is imposed upon 
the phase- squeezed beam using an electro-optic modulator. Homodyne detection is performed 
on this phase-modulated beam with an overall efficiency of rj = 0.85. The homodyne current 
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goes to the optimal feedback filter (another low-pass filter with a cutoff frequency \/2n (26)). 
Its output, (pf(t), is then shifted by ir/2 and applied on the phase of the LO beam with another 
electro-optic modulator. 

We record (f(t), I(t) and (p/(t) by an oscilloscope with a sampling rate of 100 MHz. Figure 
2 shows a typical segment of the recorded signals, plus the smoothed estimate (p s (t). The 
parameters here are k = (1.9 ± 0.1) x 10 4 rad/s, A = (5.9 ± 0.5) x 10 4 rad/s, \a\ 2 = (1.00 ± 
0.06) x 10 6 s _1 , squeezing -3.1±0.1 dB (r m = 0.36 ± 0.01) and anti-squeezing 5.1±0.1 dB 
(r p = 0.59 ± 0.01), from a pump beam power of 80 mW. Note that k and A are fixed through 
this paper. The current I{t) has zero mean because the feedback loop is designed to operate 
the homodyne measurement at this point of maximum sensitivity (Fig. IB). While the filtered 
estimate ipf(t) has a visible delay due to its causal nature, the smoothed estimate <p s (t) does 
not, and the signal phase <p(t) is reliably tracked. 

To investigate the squeezing-enhancement, we perform phase tracking with a fixed square 
amplitude \a\ 2 = (1.00 ± 0.06) x 10 6 s" 1 but with varying squeezing levels arising from OPO 
pump beam powers of 0, 30, 80 and 1 80 mW. Independently of the phase estimation, squeezing 
and anti-squeezing levels were measured for each pump beam power (26). The red crosses in 
Fig. 3 show the MSEs of the smoothed estimates a 2 as a function of the squeezing level. The 
MSE was calculated from 2 ms of data (2 x 10 5 samples). Repeating this 15 times gave the 
average MSE and its standard deviation. 

Figure 3 shows three key results. First, the squeezing-enhancement is verified: the MSEs are 
reduced below the CSL (i) by using phase-squeezed beams. Second, the experimental results are 
in good agreement with the prediction (ii), and in disagreement with the theory curve (iii) which 
is based on approximating the homodyne output current I{t) to only first order in [<p(t) — <ff{t)] 
so that _R sq = e~ 2r ' m . Third, at the higher squeezing level the MSE is saturated, indicating the 
existence of an optimal squeezing level. Even in the theoretical curve (iv) for pure squeezed 



5 



beams and zero loss, the MSE has a minimum. This curve corresponds to the fundamental limit 
imposed by Heisenberg's uncertainty principle for the phase and amplitude quadatures, namely 
e~ 2r ' p x e 2rm > 1. Although more squeezing decreases the e~ 2rm term in _R sq , it increases the 
e 2Tp term due to the tracking being imperfect, which itself is a consequence of the noise in 
the photocurrent 0. This defines (self-consistently) the optimal degree of squeezing, which 
depends upon the parameters \a\, k and A (26). 

Experimentally, we varied the amplitude |a|, while fixing the pump beam power to 80 
mW, giving squeezing and anti-squeezing levels of — 3.2±0.2 dB and 4.9±0.3 dB, respectively. 
Theoretically, the optimal squeezing level increases with \a\, and so too does the squeezing - 
enhancement, without limit. However, for our experimental conditions (10 6 s _1 < \a\ 2 < 
10 7 s _1 ) the effect of keeping the squeezing fixed is minor (less than 3% difference to a 2 ). Fig- 
ure 4 shows the dependence of the MSE a 2 on |a|. The theoretical curves show good agreement 
with experiments. Over the whole amplitude range, the estimates with the squeezed beams sur- 
pass what is possible with coherent states, with a 2 , averaged over the four different amplitudes, 
being (15 ± 4)% below the CSL. The conclusion is essentially unaltered if one calculates the 
CSL not in terms of \a\ 2 , but in terms of the effective photon flux ftf c s, which equals |ct| 2 plus 
the extra photons due to the squeezed vacuum fluctuations in the relevant spectral range (26). 

We have tracked the phase of a squeezed optical field that varies stochastically in time over 
a significant angular range. Our use of Kalman filtering in real-time adaptive measurements 
of nonclassical systems could be applied also in solid-state and nanomechanical devices. Opti- 
mizing both the degree of squeezing and its bandwidth according to the experimental conditions 
would allow a completely rigorous treatment of photon flux. Lower losses and more squeezing 
would then enable a dramatic improvement to a precision that scales differently with photon 
flux, with a 2 oc A/"~ 5 ^ 8 (20) as opposed to the a 2 oc A/" -1 / 2 in the current setup. 
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Fig. 1 Experimental setup. (A) Schematic of the experiment. (B) Homodyne output current 
versus relative phase between LO and phase-squeezed beam, and phasor diagram for a slightly 
non-optimal relative phase (as will occur in the phase tracking problem). (C) Detail of the 
experimental setup. The abbreviations are TiS: Titanium Sapphire, CW: continuous-wave, PM: 
phase modulator, SHG: second-harmonic-generation, OPO: optical parametric oscillator. 
Fig. 2 Time domain results of phase estimate. (A) Signal phase to be estimated (p(t). (B) 
Homodyne output current I(t). (C) Filtered estimate ff{t). (D) Smoothed estimate tp s (t). 
Fig. 3 Smoothed MSE a 2 versus squeezing level. Red crosses represent experimental data. 
Trace (i) is the coherent- state limit which is reachable with a coherent beam only if we have 
unit-detection efficiency rj — 1. Trace (ii) is the theoretical curve from Eq.[3l Trace (iii) is the 
theoretical curve based on a approximating the homodyne output current I(t) to only first order 
in [<p(t) — <ff{t)] so ma t -^sq = e~ 2rm . Trace (iv) is the theoretical curve from Eq. [3] for pure 
squeezed beams (i.e. without loss). 

Fig. 4 Dependence of the smoothed MSE a 2 on the amplitude squared \a\ 2 . Blue and red crosses 
are experimental data for coherent and squeezed beams, respectively. Trace (i) is the coherent- 
state limit. Trace (ii) is the theoretical curve for coherent beams with the experimental setup 
(i.e. including inefficiency). Trace (iii) is the theoretical curve for squeezed beams, including 
inefficiency. Trace (iv) is the theoretical curve for pure squeezed beams and 100% efficiency, 
with the squeezing level optimized for each \a\ 2 . 
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Supplementary Materials for 
Quantum-enhanced optical phase tracking 

SI Introduction 

Specific details of a more technical nature omitted in the main text are briefly discussed here. 
Additionally, the following sections cover a more technical discussion the experiment and de- 
tails of the mathematics supporting the arguments in the main text. 

Optical phase estimation has been extensively investigated because of its great importance 
in both science and engineering. As stated in the main text, phase estimation can be divided 
into phase sensing, where the phase is known to always lie within some small interval, and 
unconstrained or general phase estimation, where it is not so constrained (e.g. (11)). In phase 
sensing, if (as is typically the case) the field has a large coherent amplitude, the problem can 
be linearized in terms of the unknown phase-rotation (7). This means that the unknown phase- 
rotation can be treated as an unknown quadrature-displacement, considerably simplifying the 
problem. The problem that we address in our paper is one of general phase estimation, where 
the phase cannot be approximated by a fixed quadrature (see Sec. S3). Moreover we consider a 
phase which varies continuously and stochastically in time. 

To attack this problem we use an adaptive estimation technique, and sophisticated data pro- 
cessing. We estimate the phase at any time with a mean square error (15 ± 4)% below the 
coherent- state limit (CSL) — the best precision that can be achieved with classical (i.e. co- 
herent) states of light. This is possible because we use nonclassical phase-squeezed light, so 
called because the phase uncertainty is squeezed below the CSL (7-9, 19, 20). Because of 
the Heisenberg uncertainty principle, phase squeezing necessarily results in anti-squeezing in 
the amplitude. This anti-squeezing plays a key role in our experiment, causing there to be an 
optimal degree of squeezing, unlike in gravitational wave detection where more squeezing is 
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better (7-9). Our observation of this optimal squeezing is a new phenomenon in quantum phase 
estimation. 

S2 Details of the experiment 

Figure SI shows the experimental setup in more detail than shown in Fig. 1C in the main text. 
Specifically we show an optical mode cleaner, acousto-optic modulators (AOMs) and a radio 
frequency (RF) source, all of which were omitted from Fig. 1C in the main text in the interests 
of clarity. 

The signal (p(t), which is the waveform to be estimated, is imposed upon a weak phase- 
squeezed state using an electro-optic modulator. Adaptive homodyne detection of the phase- 
modulated beam is used to form the phase estimate. 

In this section we describe the generation of the weak phase- squeezed state, and then de- 
scribe its characterization using static homodyne detection. We then describe the generation of 
the signal and the phase estimation procedure using adaptive homodyne detection. 

S2.1 Generation and measurement of phase-squeezed states 

We use an 860 nm continuous-wave (CW) Titanium Sapphire laser in the experiment. A fraction 
of the laser beam is used as the local oscillator (LO) for homodyne detection, with a spatial- 
mode cleaner used to improve mode-matching of the LO to the phase-squeezed state. The 
phase- squeezed state is generated by a sub-threshold optical parametric oscillator (OPO) using 
a periodically poled KTiOP04 crystal (27). The OPO is pumped by the 430 nm output of the 
second-harmonic-generator (SHG) cavity. The AOMs are used to generate a pair of optical 
sidebands at ±5 MHz (21), which are equivalent to a weak coherent state. The optical side- 
bands, which are within the 13 MHz half-width-half-maximum bandwidth of the OPO, are then 
injected into the OPO. 
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Figure SI: Experimental setup. The abbreviations are TiS: Titanium Sapphire, CW: 
continuous-wave, AOM: acousto-optic modulator, PM: phase modulator, SHG: second- 
harmonic-generation, OPO: optical parametric oscillator, RF source: radio frequency (5 MHz) 
source. 

The relative phase between squeezing and the 5 MHz sidebands must be locked to ensure 
that a phase-squeezed beam is generated. For that purpose, we tap off 1% of the OPO output 
beam to be used as a signal to lock the relative phase between the squeezing and the sidebands. 
Since the OPO output power is extremely low (on the order of pW), we use another homodyne 
detector to measure the tapped beam (not shown in Fig. SI). The method is similar to that used 
in Ref. (27). 

In addition to the deliberately-imposed phase modulation signal (p(t), the experiment is 
subject to low-frequency phase disturbances from the environment. We use a low-frequency, 
low-gain feedback loop to suppress the environmental disturbances. When we use static homo- 
dyne detection to characterize the squeezing (anti-squeezing) level of the phase-squeezed states 
before performing the phase estimation experiment, we use only this classical feedback loop to 
lock the relative phase between the LO and the phase- squeezed beams. 

With a LO beam power of 2 mW, the shot noise level of the homodyne detector is 10.4 dB 
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higher than the circuit noise level (that is, the ratio of the measured shot noise to circuit noise 
is S = 11.0 ). The homodyne visibility is £ = 0.988, the optical transmission is p = 0.97, and 
the quantum efficiency of each of the photo-diodes is £ = 0.99. The overall efficiency of the 
homodyne detection is thus calculated as rj = £ 2 p((S — 1)/S = 0.85 (28). 

Figure S2 shows static homodyne measurements of the phase-squeezed state when the 
power of the OPO pump beam is 180 mW. Fig. S2 shows the noise spectrum as measured 
with a radio frequency spectrum analyzer in the region around the 5 MHz sidebands. As we 
have mentioned previously, these are static homodyne measurements of the phase- squeezed 
state. So for Fig. S2 we do not apply phase modulation signal tp(t) to the phase- squeezed beam, 
and do not use the phase-tracking loop. The residual peak at around 5 MHz in trace (ii) comes 
from fluctuations of the classical locking loop. Compared to our signal bandwidth (\/2n = 9.4 
kHz), the bandwidth of the squeezing is broad enough that we may assume that the squeezing 
is broadband over the frequency range of interest. 

In our experiment, we use pump beam powers of 0, 30, 80 and 180 mW. Squeezing and 
anti-squeezing levels around 5 MHz are shown in Fig. S3. The solid line is a theoretical curve 
fitted to the experimental results with the overall loss l sq — 1 — 7/(1 — l a ) (comprising the 
homodyne detection efficiency, 77, and the loss from the OPO and phase modulator, l s ) being a 
free parameter. We define the measured squeezing and anti-squeezing levels as R_ := e~ 2Vm 
and R + := e 2r ' v . With an overall loss / sq and a pure squeezing parameter r, the measured 
squeezing and anti-squeezing levels can be written as (29), 
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In Fig. S3, we obtain the overall loss Z sq = 0.33 which mainly comes from the OPO and phase 



modulator. 
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Figure S2: Characterizing the phase squeez- 
ing. The data is recorded with a spectrum an- 
alyzer. Resolution and video bandwidths are 
10 kHz and 300 Hz, respectively. Trace (i) is 
shot noise level. Traces (ii) and (iii) are the 
squeezing and anti-squeezing levels where 
the LO phase is locked to those quadratures. 
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Figure S3: Anti-squeezing level versus 
squeezing level. Solid curve is fitted to the 
experimental data with loss being a free pa- 
rameter. The overall loss l sq is obtained as 
0.33. 



S2.2 Phase Tracking 

The signal waveform to be estimated is generated by a random signal source and a low-pass 
filter with a cutoff frequency of X/2-k = 9.4 kHz. The signal source is a digital signal genera- 
tor (AFG3021, Tektronix), which generates white Gaussian noise (3dB bandwidth is 25 MHz). 
The signal ip(t) is imposed upon the phase-squeezed beam using an electro-optic modulator 
(LM0202, Linos). We perform homodyne detection on this phase-modulated beam. After 
homodyne detection of the phase-modulated beam, the homodyne current is demodulated at 
5 MHz and then goes to the tracking filter (another low-pass filter with a cutoff frequency 
A/27T according to the Kalman filter theory). The output of the tracking filter gives the filtered 
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estimate <Pf(t), which is then shifted by ir/2 and applied to the phase of the LO beam with 
another electro-optic modulator (waveguide modulator, EOSPACE ). The gain of the phase- 
lock loop is optimized to give the maximum sensitivity. As explained previously, an additional 
low-frequency, low gain feedback loop is used to isolate the experiment from environmental 
disturbances. 

We record (p(t), y?/(t) and the demodulated homodyne current I(t) using an oscilloscope 
with a sampling rate of 100 MHz. We calculate the mean square error from 2 ms of data, 
comprising 2 x 10 5 data points. To calculate the average and the standard deviation of the 
mean square error, we repeat the measurement 15 times. Note that each the correlation time 
of the noise signal cp(t) is approximately A -1 , and so 2 x 10 5 data points have effectively 
2 ms/A^ 1 ~ 10 2 independent data points. 

S3 Phase estimation with broadband squeezing 

In this section we describe phase estimation with broadband squeezing. As in the main text, the 
phase variation ip(t) obeys the following, 



where dV(s) is a Wiener increment which satisfies (dV(s)dV(s')) = S(s — s')(dt) 2 , A is the 
bandwidth of the phase noise if, and k/2\ is the mean square variation of <p. Note that in our 
experiment this variation is much larger than mean square error with which we can estimate 
cp. The phase-modulated beam is measured by homodyne detection. The LO phase $(£) is 
adaptively controlled to be $(£) := ipf(t) + tt/2, where (pAt) ~ ip(t) is our filtered estimate. In 
our experiment the squeezing bandwidth is many MHz, which is much larger than the bandwidth 
A of the phase noise. Hence we can approximate the squeezing noise as white. Using the 




(6) 



20 



normalization of (19), the homodyne output current is 

I{t)dt = 2 \a\ sin [<p(t) - <p f (t)] dt + J ' R sq {t)dW{t), (7) 
R sq (t) = sin 2 [<p(t) - <p f (t)]e 2r ' + cos 2 [<p(t) - <p f(t)]e~ 2r ™, (8) 

where |a| 2 is the photon flux of the coherent component of the beam (i.e. photons per unit time), 
dW(t) is a Wiener increment representing quantum white noise, satisfying (dW(s)dW(s')) = 
S(s — s')(dt) 2 , and the amplitude of this noise R sq is determined by the squeezing (r m > 0) 
and anti- squeezing (r p > 0) parameters defined earlier. The squeezing factor R sq (t) is time- 
dependent because of tracking error p(t) — (pf(t), the unknown and varying misalignment be- 
tween the squeezing ellipse and the local oscillator phase. 

We denote the tracking error by A/(t) = ip{t) — ipf(t), and its stationary ensemble average 
by cr 2 = (A 2 .) . Because Of 1 in our experiment we can approximate I{t)dt by a 2nd order 
expansion in Af(t). This gives 

I(t)dt ~ 2|o;| A f (t)dt+ y/R^dW(t), (9) 
it! sq (t) ~ A5(t)e^ + [1-A 2 f (t)]e- 2 ^. (10) 

The approximations given here help to illustrate the difference between phase estimation 
and phase sensing. In phase sensing, the variation in the system phase (p(i) is sufficiently 
small that I{t)dt can be linearised in ip(t). In phase estimation, this is not the case and a 
better approximation is needed. Here the filtered estimate <ff(t) is sufficiently good that the 
sine function can be linearised in A/(i). (This is why feedback is vital in phase estimation.) 
Furthermore, in the case where there is squeezing, there is a difference in the treatment of 
the squeezing term R sq (t). In phase sensing, the variation in the system phase is sufficiently 
small that one may simply take R sq (t) ~ e~ 2rm . One can simply treat the phase-rotation as 
a displacement of the squeezed quadrature, and not consider the anti- squeezing. In contrast, 
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for phase estimation (even with feedback) the variation of the system phase is sufficiently large 
that one needs a better approximation for R sq (t). Because the filtered estimate is not perfectly 
accurate, there is a significant contribution to the error from the anti-squeezed quadrature. 

A further approximation is to replace A*(i) by its stationary average a 2 , which is indepen- 
dent of time. (This approximation is justified below by comparing the time scale of its variation 
to that of the subsequent filter.) The normalized homodyne output current can then be written 
as 

I(t)dt ~ 2\a\A f (t)dt+ JR^dWit), (11) 
R sq = a 2 f e 2r - + (l-a 2 f )e- 2r -, (12) 

where R sq is an effective (time-independent) squeezing factor which takes into account the 
tracking error in an average sense. A measurement better than the coherent- state limit is ex- 
pected to be possible when this _R sq is less than 1. 

For an equation of the form of Eq. QT1 it is known that the optimal causal estimate (as 
required for tracking) is provided by the Kalman filter (22): 

<p f (t)=T f e^^ds, (13) 

J-oo 2 I«I 

where T is the Kalman gain. Before going further, we return to the approximation in Eq. [T2l 
where we replace the time- varying A^(£), with its average a 2 . The characteristic time of the 
Kalman filter is A -1 ; that is, the homodyne output current I(t) is, roughly speaking, averaged 
over the period of A -1 . Thus if the variation of A^(£) is rapid compared to A -1 , it is justified 
to replace the time- varying A 2 (i) by the average a 2 . Figure S4 shows a typical experimental 
trace of the temporal variation of the squared filtered estimate error A 2 Jt) = [ip(t) — pf(t)} 2 , 
and also its auto-correlation function. The experimental conditions are the same as Fig. 2 in the 
main text. The variation of the squared filtered estimate error has a correlation time around 2/xs. 
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The Kalman filter characteristic time is calculated as A~ x ~ 17/is which is much larger than the 
correlation time of the squared filtered estimate error. Therefore our approximation in Eq. [12] 
turns out to be appropriate. 




T (US) 

Figure S4: (A) Squared error of the filtered estimate Ay under the same conditions as Fig. 2 
in the main text. The purple trace is A|, and the blue trace is that after the Kalman filter. (B) 
Normalized auto-correlation ([A/(t)] 2 , [Af(t + r)] 2 ) as measured in the experiment (purple 
curve), and theoretical curve for the Kalman filter of white noise (blue curve). 



This separation of time scales is a consequence of the fact that our experiment was per- 
formed in the limit <C A -1 (for the traces in Fig. S4, T^ 1 « 4/xs). This can be seen by 
considering the differential equation for Af(t). First, combining Eqs. [TT1 and [T3l gives 



Combining this with Eq. [6] gives 



t 

I e 

oo 



-A(i-s) 



A f (s)ds 



2 \a\ 



dW{s) 



(14) 



dA f = d<p(t) - d<p f (t) 



-(A + T)A f (t)dt + y/RdV{t) - T yrf* dW(t). 



(15) 
(16) 
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Thus, under the above approximation, A f has a correlation time of (A + r) _1 , and A 2 a correla- 
tion time of half that. For T ^> A, the correlation time for A 2 is much less than that of the filter, 
and hence on the time scale of the filter A 2 may be replaced by its ensemble average. More 
fundamentally, this difference in time scales comes because we have conducted our experiment 



in the regime where e = J X 2 R sq / (4 | or | 2 k) is small, as in this regime r ~ A/e (see below). 

Even if we do not replace R sq (flOl) by _R sq (fl~2l . we find the same result for er?. This follows 
from the differential equation for (A 2 ) using Eq.[lOl namely 

T 2 (A 2 e 2rp + (1 — A 2 )e~ 2rm ) 
d(A 2 f ) = -2(\ + T)(A})dt + Kdt ^ -±- 2 f - '-dt. (17) 

This gives a linear equation for a 2 = (A 2 ) ss which is identical to the equation which would be 
obtained from Eq.[T6]in which R sq contains a 2 already. This means that our theory models the 
experiment for any value of e (provided a 2 remains small), but our argument for the optimality 
of the Kalman filter, and the corresponding smoothed estimator, holds only for e small. This is 
reasonable for our experiment, in which e < 0.2. 

Returning to Eq.[[4l it is obvious that phase estimation using the phase- squeezed beam with 
the amplitude \a \ and squeezing factor R sq is mathematically equivalent to that using a coherent 
beam with the amplitude \a\/^R sq . (Note that this equivalence holds in terms of the estimate 
precision, not in terms of photon number.) Therefore we can obtain the Kalman gain T and the 
mean square error a 2 by formally replacing |a| in Ref. (22) with |a|/\/ -R S q> giving 



T = -A + ^X 2 + 4K\a\ z /R sq , (18) 

In the limit e 1, the optimal gain is T w A/e, and a 2 w a/ nR sq / (2\a\). All of these 
equations are still implicit because -R sq (fl~2l) is a function of a 2 . Solving the implicit equation 
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for a 2 f gives explicitly, and exactly (i.e. without assuming e<l) 



-A - 



k (e 2rp — e 2r " 



2 e -2r„ 

1 



k (e 2r p — e _2r " 



-2r„ 



2A 



16|a| 2 K 

— 2r ' 



(20) 



| a | 2 + 4A(e 2r p — e" 



•2r„ 



-2Ae- 2rm + K (e 2rp - e" 2r " 



+ y [2Ae _2rm + n{e 2r p — e~ 2Tm )} + 16|a| 2 Ke" 



2r n 



(21) 



As in Ref. (27), we can apply the smoothing technique of Ref. (22) for further improvement. 
The smoothed estimate <p s (t) and mean square error a 2 = (Up(t) — ip s (t)} 2 } are 



<p s (t) = (2A + T) / e-^^iffis)^ 



H 



(22) 
(23) 



2y / 4K|a| 2 /# sq + A 2 

These formulae are obtained by formally replacing \a \ in Ref. (22) with \ a\/ ^/ R sq - In the limit 
e 1, we obtain the simple relation a 2 w cr 2 /2 as in Ref. (2i). 

S4 Optimum squeezing level for phase tracking 

In this section, we discuss the optimum squeezing level for the phase tracking, as quantified by 
the mean square error in the smoothed estimate. 

From Eq.[T2]and[23l it is obvious that highly squeezed states would not improve the phase 
tracking if there were too much anti- squeezing noise. Quantum mechanics does not allow us 
to squeeze a physical quantity without anti- squeezing the canonically conjugate counterpart. 
Hence there must be an optimum squeezing level for the phase estimation problem. Figure S5 
shows the smoothed mean square error a 2 as a function of both squeezing and anti-squeezing 
levels, with other parameters (a, k, A) fixed. In this figure the lower-right half is forbidden due 
to the uncertainty principle. Even for a pure squeezed state (at the boundary of the forbidden 
region), the minimum mean square error is obtained for an optimum squeezing level of about 
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7 dB. In practice, with finite loss, the squeezed state is not pure, and it is not possible to attain 
the boundary with the forbidden region. The minimum anti-squeezing for a given squeezing 
will follow a curve away from the boundary, according to Eq. [5] The curve for the loss in our 
experiment is shown in Figure S5. The qualitative behavior is the same as in the pure state case. 
There is an optimal squeezing level, and increasing the squeezing beyond this will increase the 
mean square error. 




Squeezing level (dB) 

Figure S5: Mean square error of the smoothed estimate, as a function of both squeezing and 
anti-squeezing levels, calculated from Eq. [23] Other parameters (a, k, A) are fixed and the 
same as Fig. 3 in the main text. The black dashed line and red crosses are the predicted and 
experimentally observed squeezing and anti-squeezing levels in our setup (see Fig. S3). CSL 
stands for the coherent- state limit. 

The optimum squeezing level varies with parameters a, k and A. From Eq. [T21 it would 
be obvious that the optimum squeezing level is higher if the mean square error of the filtered 
estimate a'j is smaller. This is because the anti-squeezing noise has less effect if cr| is small. 
The mean square error aj would be reduced by increasing a, or A, or by reducing k. Note that 
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the mean square variation of signal <p is k/2\, and the estimate error aj is smaller with smaller 
signal variation (smaller k or larger A). Therefore, the optimum squeezing level is higher if a 
or A is larger or k is smaller (see Fig. S6). 




io 4 io 6 io 8 io 4 io 5 io 6 io 3 io 4 io 5 
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Figure S6: Optimum squeezing levels as a function of (A) \a\ 2 , (B) A and (C) k. The optimum 
squeezing levels are calculated by minimizing a 2 of Eq. |23j The parameters other than the 
varying parameter are fixed and the same as Fig. 3 in the main text (e.g., in (A), k and A 
are fixed and the same as Fig. 3 in the main text.). Traces (i) are for squeezed beams in our 
experimental setup (i.e. including loss). Traces (ii) are for pure squeezed beams. 



S5 Squeezing bandwidth and photon flux 

In this section, we derive the photon flux due to squeezing, as this is required (see Sec. S7) 
for a complete accounting of the resources used in phase estimation. As we shall show shortly, 
infinite-bandwidth squeezing leads to an infinite photon flux. Therefore, finite-bandwidth squeez- 
ing must be considered explicitly. 

Let us consider a squeezing spectrum R_(Q) and anti-squeezing spectrum R + (£l) of the 
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standard form (27) 




(24) 



x 



R+-R-- 2y/(l -R_)(R+ - 1) 
i? + + i?_ - 2 



(25) 



where R± are squeezing and anti- squeezing levels at the center frequency (Q = 0). In the 
case of the OPO, 2Af2 and x correspond to the cavity's decay rate and the normalized pump 
amplitude respectively. The difference between the bandwidths for squeezing (1 + x) Af2 and 
anti- squeezing (1 — x) Af2 ensures that R + (VL)R_(VL) = 1 for all Q when the squeezed state is 
pure. 

The number of photons due to squeezing in a frequency mode f2 is given as (29), 



The photon flux per unit time J\f sq due to the squeezing is given by integrating over (30), 



Note that, as mentioned previously, the photon flux diverges for infinite-bandwidth squeezing. 

S6 Phase estimation with finite-bandwidth squeezing 

In this section, we derive mean square error of the phase estimate when the squeezing bandwidth 
is finite. This will be used in Sec. S7 to calculate the (small) effect that restricting the bandwidth 
of our squeezing would have had. 

For finite-bandwidth squeezing the equations QT] and [121 for the homodyne output current 
must be modified to 



n(n) = ^[i2 + (n) + i2_(n)]-i 



(26) 



1 f°° 1 Afi 

Af sq = / n(Q)dQ = - [{R+ - 1)(1 - x) + {R- -!)(! + x)] (27) 




I(t)dt 



2\a\ [tp{t) -<p f (t)]dt + dY(t) 



(28) 




(29) 
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where X±(s — s') are auto-correlation functions for the squeezed (— ) and anti-squeezed (+) 
quadratures, and the tilde is used to distinguish particular finite-bandwidth properties from the 
previously calculated infinite-bandwidth quantities (with no tilde). X±(s — s') are calculated 
by the inverse Fourier transform of the squeezing and anti- squeezing spectrum R±(Q), 

X ± {s-s') = — / R ± {Q)e inis - s,) dQ 

= 6{s _ s >) + {R± _ 1) ( 1 I ^) A ^V (i T .)AO |W| (30) 

From Eqs. [T3]and[2H and assuming a filter of the same form as above, the estimator used 
for tracking is 



J— oo ^1^1 J— oo I 



dY(s) 
2 bl 



(3D 



where the optimal Kalman gain V can be found numerically. In order to solve this equation, we 
transform it into differential form as 



dip f {t) = -\tp f (t)dt + rU(p(t)-0 f (t)]dt 



dY(t) \ 



2\a\ 



-{x + r)0 f {t)dt + f 



<p(t)dt + 



dY(t) 
2\a\ 



(32) 



Then we transform it into integral form again, 



0f(t) 



-(A+f)(t-s) 



ip(s)ds 



dY(s) 
2\a\ 



(33) 



From this, we can see that the noise term dY(t) is low-pass filtered with a cutoff (angular) 
frequency (A + T). Thus we could consider that the bandwidth of the estimator is (A + T). 
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We calculate the mean square error aj as, 
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(34) 
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Here (<£>(s)</?(s')) is calculated from Eq.[6]as, 



(35) 



Now a? can be given from Eqs. [2<3[30l[34]and[35l 



a 2 
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i r 2 



2(A + T) 8|a| 2 r + A 

h 4 



x |l + a 2 f (i? + - 1 
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h+ + l 
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(36) 
(37) 



A + r 

The effect of finite-bandwidth squeezing is characterized by h±, which is the ratio of the squeez- 
ing (anti-squeezing) bandwidth (1 + x)Ail to the estimator bandwidth (A + T). In the limit 
Af2 — > oo we find from Eq. [36] that cr 2 obtains its minimum value of cr 2 when T = T, as 
expected. 

The optimal smoothed estimate with finite-bandwidth squeezing (p s {t) is similarly found to 

be 

(p s {t) = (2A + f) e'^ f ^(p f {s)ds, (38) 
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and the mean square error of the smoothed estimate is 

_ 2 k(2\ 2 + f 2 + 2Af) f 2 (2A + f) 2 
4(A + f) 3 + 16|a| 2 (A + f) 3 

Similarly to that for the filtered estimate, the effect of finite-bandwidth squeezing is character- 
ized by h±. Again, in the limit Af2 — > oo we find from Eq. [39] that a 2 obtains its minimum 
value of a 2 when T — V. 

S7 Effective squeezing bandwidth and photon flux 

In our experiment, we use narrow optical sidebands as our signal. The frequency range of 
squeezing which is exploited, in the vicinity of the sideband frequency, is thus much narrower 
than the OPO's bandwidth. In other words, most of the photon flux produced by the squeezing 
plays no part in the current experiment. Counting the entire photon flux of the beam would 
grossly misrepresent the photon flux resource used by our phase tracking algorithm. To reason- 
ably represent the photons in the squeezing (in addition to those in the coherent field) we must 
consider the photons in a relatively narrow band in the vicinity of the sideband frequency. That 
is, we must define the effective squeezing (and anti- squeezing) bandwidths. 

Once again we model the spectrum of the squeezing and anti- squeezing by Eq. [2H with 
£1 = denoting the sideband frequency. As before, R± denotes the squeezing (anti-squeezing) 
levels at the sideband frequency, but now we replace Af2 by an effective bandwidth Af2 eff to 
model narrow-band squeezing. The effective bandwidth is the narrowest bandwidth of squeez- 
ing that could have been used without substantially affecting the experimental results obtained. 
A reasonable guess would be to set the effective bandwidth somewhat wider than the estimator 
bandwidth. This is because the squeezing within the estimator bandwidth can effectively en- 
hance the estimate, whereas the squeezing outside of the estimator bandwidth is mostly wasted. 
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Indeed, it turns out that defining Af] c fj = 2(A + T) is appropriate in the sense that a 2 is little 
different from a 2 in the regime of our experiment. 

From now on, we restrict to the experimental parameters as shown in Fig. 4 of the main text. 
The squared coherent amplitude |ct| 2 ranges from 10 6 to 10 7 s _1 , and the other parameters are 
fixed as k = 1.9 x 10 4 rad/s, A = 5.9 x 10 4 rad/s, iL_ = 0.479, and R+ = 3.09 (see the main 
text). The normalized pump amplitude x is calculated as x = 0.33 from Eq. [251 In Fig. S7, 
we plot the mean square error of the phase estimate as a function of \a\ 2 (remember that T is a 
function of \a\ 2 in Eq. [20l and hence AQ e s is too). The green line indicates the mean square 
error a 2 for broadband squeezing (Af2 eff — > oo). The red line is the mean square error a 2 with 
Af2 c fj = 2(A + T). The difference between a 2 and a 2 is less than 3%, which is smaller than the 
experimental error bars which are ~ 10%. 

By setting Af] cff = 2(A + T), we can now calculate photon flux of the squeezing J\f^ from 
Eq.[27l and obtain the total effective photon flux M cS = \a\ 2 + Af^. The contribution of 
to the total photon flux is at most 7% under our experimental conditions. Figure S8 shows the 
mean square error versus the effective photon flux J\f e g under the same conditions as Fig. 4 in 
the main text. We also replot the experimental data from that figure. To make an absolutely 
fair comparison, we correct the experimental mean square errors for broadband squeezing to 
the expected (i.e. slightly larger) values for the effective bandwidth. The key point is that even 
taking into account the photon flux due to squeezing in the effective bandwidth around the 
sidebands, the accuracy of our phase estimates surpasses the ideal (lossless) coherent- state limit 
as a function of total photon flux over the whole range of the experiment. 
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Figure S7: Theoretical curves for the mean 
square error in the smoothed estimate ver- 
sus the photon flux | or | 2 due to the coherent 
amplitude. Experimental conditions (k, A, 
and R±) are same as Fig. 4 in the main text. 
The green line is erf for broadband squeez- 
ing (as in the experiment). The red line is a 2 s 
for finite-bandwidth squeezing with Af2 cff = 
2(A + r), where T is the feedback gain. 
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Figure S8: Mean square error of the 
smoothed estimate versus total effective pho- 
ton flux Atff. Parameters are as in Fig. S7. 
Blue and red crosses are experimental data 
for coherent and squeezed beams respec- 
tively. The latter are corrected to the ex- 
pected values for squeezing with bandwidth 
Af2 cff . Trace (i) is the theoretical coherent- 
state limit. Trace (ii) is the predicted (i.e. 
with experimental imperfections) curve for 
coherent beams. Trace (iii) is the predicted 
curve for finite-bandwidth squeezed beams 
including all experimental imperfections. 
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